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Abstract 

In this paper we analyze entanglement classification of extended Greenberger-Horne-Zeilinger- 
symmetric states p ES , which is parametrized by four real parameters x, yi, y2 and y%. The 
condition for separable states of p E is analytically derived. The higher classes such as bi-separable, 
W, and Greenberger-Horne-Zeilinger classes are roughly classified by making use of the class- 
specific optimal witnesses or map from the extended Greenberger-Horne-Zeilinger symmetry to the 
Greenberger-Horne-Zeilinger symmetry. Prom this analysis we guess that the entanglement classes 
of p ES are not dependent on yj (j = 1, 2, 3) individually, but dependent on yi + y2 + y3 collectively. 
The difficulty arising in extension of analysis with Greenberger-Horne-Zeilinger symmetry to the 
higher-qubit system is discussed. 
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I. INTRODUCTION 

Entanglement [l| is an important physical resource in the context of quantum information 
theories^]. As shown for last two decaces it plays a crucial role in quantum teleportationjsj, 
superdense coding^], quantum cloning 5], quantum cryptography [6]. It is also quantum en- 
tanglement, which makes the quantum computer outperform the classical one|7(]. Therefore, 
it is greatly important task to understand what kind and how much entanglement a given 
quantum state has. 

For multipartite quantum states there are several types of entanglement. Each type is in 
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general categorized by stochastic local operations and classical communication (SLOCC) 
Thus, these types of entanglement is often called SLOCC-equivalence classes. For example, 
for three-qubit pure states(9j] there are six SLOCC-equivalence classes such as separable, 
three bi-separable (A — BC, B — AC, C — AB), W and Greenberger-Horne-Zeilinger (GHZ) 
classes. Among them genuine tripartite entanglement arises in W and GHZ classes. The 
representative states of these classes are 

|GHZ) = i=[|000) + |lll)] (1.1) 
|W) = i=[|001) + |010) + |100)]. 

One of the most remarkable fact in this classification is that the set of W states forms 
measure zero in the whole three-qubit pure states. 

This classification can be extended to three-qubit mixed statesflO]. Following Ref.fiol] 
the whole three-qubit mixed states are classified as separable (S), bi-separable (B), W and 
GHZ classes. These classes satisfy a linear hierarchy S C B C W C GHZ. One remarkable 
fact, which was proved in this reference, is that the W-class 1 is not of measure zero among 
all mixed-states. 

Although SLOCC classes for three-qubit system are well-known, we still do not know how 
the entanglement is classified in multi-qubit system except four-qubit pure states, where 



there are nine SLOCC classes 12]. Furthermore, still it is very difficult problem to find a 



SLOCC class of a given three-qubit mixed state 2 except few rare cases. Thus, it is important 



1 As Ref. we will use the names of the SLOCC classes in an exclusive sense throughout this paper. 

2 For three-qubit pure states it is possible to find the SLOCC classes by computing the concurrence [13j of 



the reduced states and three-tangle[14( for the given states 



2 



task to develop a method, which enables us to find a SLOCC class of an arbitrary three-qubit 
states. 

Recently, a significant progress is made in this issue. In Ref.[15| a complete SLOCC 
classification for the set of the GHZ-symmetric states was reported (see Fig. 1). According 
to this complete classification the ratio of number of S, B, W, and GHZ states in the whole 
set of the GHZ-symmetric states is 1 : 1 : 1.076 : 0.924. Thus, W class is not of measure 
zero in this restricted set of the three-qubit states. Using this classification the three-tangle 
r for the arbitrary GHZ-symmetric states is explored in Ref.[ii|. Moreover, this complete 
classification is used to construct the class-specific optimal witnesses for the three-qubit 
entanglement . 

The purpose of this paper is to explore a possible extension of Ref.[l5j to treat more 
three-qubit mixed states. For this purpose we enlarge the symmetry group to, so-called, 
the extended GHZ symmetry group. The whole set of quantum states invariant under the 
extended GHZ symmetry group is parametrized by four real parameters (x, y±, y2, Jte). The 
complete classification for S states is analytically derived. However, the classification for B, 
W, and GHZ states is incomplete. Rough classification for B, W, and GHZ states is explored 
by making use of the elass-specinc optima. witnessesQ or GHZ ~yQ, respeetive.y. 

The paper is organized as follows. In next section we review Ref. [15J briefly. In section 
III we discuss on the extended GHZ symmetry. It is found that the set of extended GHZ- 
symmetric states is parametrized by four real parameters. In section IV we derive a condition 
for the separable region in the four-dimensional parameter space by applying a Lagrange 
multiplier method. In section V we perform a entanglement classification of the extended 
GHZ-symmetric states roughly by making use of the class-specific optimal witnesses or map 
from extended GHZ symmetry to GHZ symmetry. In section VI a conclusion is given. In 
particular, we discuss on the difficulty arising when we extend the analysis with the GHZ 
symmetry to the higher-qubit systems in this section. 



II. CLASSIFICATION OF GHZ-SYMMETRIC STATES 

In this section we review Ref. [jj]] briefly. The GHZ-symmetric states are the three-qubit 
states which are invariant under the following transformations: (i) qubit permutations, (ii) 
simultaneous three-qubit flips (i.e., application of a x £g> a x <E> cr x ), (iii) qubit rotations about 
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the z-axis of the form 



U((j>i,fa) = e i01fTz ® e i<fe<Tz ® e -^ 1+<fe K (2.1) 

It is straightforward to show that the general form of the GHZ-symmetric states p s is 
parametrized by two real parameters x and y as 

p 5 (x )? /) = [ x+ — y + - ) |GHZ+)(GHZ + | + [ + — y+- I |GHZ_)(GHZ_| (2.2) 



+ 

where 



1 2/ 



8 2^ 



2 " 8 / ' v Tl \ 2 

|001)(001| + |010)(010| + |011)(011| + |100)(100| + |101)(101| + |110)(110| 



|GHZ±> = -^(|000>±|111». (2.3) 

The real parameters x and y are introduced such that the Euclidean metric in the (x, y) 
plane coincides with the Hilbert-Schmidt metric d(A, B) 2 = |tr(A — B)^(A — B), i.e., 

d 2 (p s ( Xl , yi),p s (x 2} y 2 )) = (x 2 - xi) 2 + (y 2 - yi ) 2 . (2.4) 

Since p s (x, y) is a quantum state, the parameters x and y are restricted as 

1 V3 2 1 , , 



4V3 ~ 4 ' y ~ v/3 4^3' 
This restriction can be easily derived by computing the eigenvalues of p s . Thus, the set 
of the GHZ-symmetric states are represented as a triangle in (x, y) plane as Fig. 1 shows. 
Each point inside the triangle corresponds to each GHZ-symmetric state. 

In order to classify the GHZ-symmetric states it is worthwhile noting that there exists a 
map from an arbitrary three-qubit pure state to the GHZ-symmetric state p s (ii), which 
is defined as 

p s (V0= / dUU\^)(^\U\ (2.6) 



where the integral is understood to cover the entire GHZ symmetry group. If, for example, 
l^) = Yli j k=o ^ijklijk), the corresponding p s {ip) is given by Eq. ( 12. 21) with 

^ = ^(^000^111 + ^000^111) V = -j= (j^oool 2 + IV'ml 2 - ^ • (2.7) 

Another fact we will use for classification is that applying GL(2, C) transformations to any 
qubit does not change the entanglement class of a multiqubit state. This fact leads from 




the invariance of the entan 
communication (SLOCC)[8 



ement class under the stochastic local operations and classical 

S|- 

In order to find the boundary of each entanglement class, therefore, we fix the y coordinate 
and derive the maximum of |x| by using Eq. (12.71) and applying the Lagrange multiplier 
method. Since mirror symmetry implies x min = —x max , it is possible to restrict ourselves to 
x > 0. This procedure yields some region in the (x, y) plane. If this region is a convex, the 
set of states corresponding to this region exhibits a same entanglement property. If it is not 
convex, the proper boundary of the class is obtained by the convex hull of this region. 

The entanglement classification derived in this way is summarized in Fig. 1. Recently, this 



classification was used to compute the three-tangle [14[ of the entire GHZ-symmetric states 
analytically 16[ . More recently, this is used to derive the class-specific optimal witnesses for 
three-qubit entanglement 111 ]. 



III. EXTENDED GHZ SYMMETRY 



In this section we will relax the condition of the GHZ symmetry to treat more large set 
of the three-qubit states. The symmetry we consider is identical with the GHZ symme- 
try without first condition, i.e., qubit permutations. We will call this the extended GHZ 
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symmetry. 

It is not difficult to show that the general form of the extended GHZ-symmetric states is 
parametrized by four real parameters x, yi, y 2 , and y 3 as 



P (x,yi,y2,y 3 ) 

i vi + m 



2/i + 2/2 — 2/3 



" 3 + x) \GHZ,){GHZ + \+ (I + Vl+V : + V3 - x 



.1 _ -2/i + 2/2 + 2/3 



|001)(001| + |110)(110| 



1 _ 2/i ~ 2/2 + 2/3 



(3.1) 

GHZ_){GHZ_\ 
|010)(010| + |101)(101| 



|011)(011| + |100)(100| 

The parameters are chosen so that the four-dimensional Euclidean metric coincides with the 
Hilbert-Schmidt metric, i.e., 



d 2 [ p ES (x,y 1 ,y 2 ,y 3 ),p I1J »(x,y 1 ,y 2 ,y 3 ) ) = {x-xf + {y x -y x f + (y 2 - y 2 f + (y 3 -y 3 f. (3.2) 



Since p ES should be a physical state, the parameters are restricted as 



, 1 1 
1 2/i + 2/2 1 - | < 2/3 < ^ - 1 2/i 



2/2 1 



n . 1 , 2/i + 2/2 + 2/3 , . , 
< — I ± x < 1. 



(3.3) 



The restriction (13.31) can be depicted pictorially. As Fig. 2(a) shows, the physically available 
value of yi {i = 1,2,3) is confined inside polyhedron in the three-dimensional (yi,y 2 ,y 3 ) 
space. As this figure exhibits, y^'s are restricted by — | < yi,y 2 ,y 3 < \. However, as Fig. 
2(b) shows, 2/i + 2/2 + 2/3 is restricted by — \ < y\ + y 2 + y 3 < | depending on x. Unlike 
the GHZ-symmetric case each point inside the triangle in Fig. 2(b) corresponds to infinite 
number of quantum states with same 2/1 + 2/2 + 2/3- 

Similarly to GHZ symmetry there exists a mapping from a set of the three-qubit pure 
states to the set of the extended GHZ-symmetric states. Let \tp) be an arbitrary three- 
qubit pure state. Then, the corresponding extended GHZ-symmetric state is given by Eq. 
( I2.6p . The only difference is a change of the symmetry group from the GHZ symmetry to 
the extended GHZ symmetry. If, for example, = Yli,jk=o ^ijkfijk) ■, the corresponding 
p E (ifi) is given by Eq.f l3.ll) with 



x 



■X (^000^111 + ^ooo^m) 



(3.4) 



2/1=2 



ooo 1 2 + |-0m I 2 + IV'oiiI 2 + l^iool 2 ) - t 



2/2 
2/3 



( I -0000 1 2 + 1-0111 1 2 + |01Ol| 2 + |0O1O| 2 ) - 
( 1 0000 1 2 + |011l| 2 + l^liol 2 + 10001 1') - 



y\ + yi + y3 




X 



FIG. 2: Pictorial representation of Eq. (|3.3p . Unlike the GHZ-symmetric case each point inside 
the triangle in (b) corresponds to infinite number of quantum states with same y\ + y 2 + 2/3. 



It is worthwhile noting a relation 



(u <g> u ® u)p ES (x, yi, y 2 , y 3 )(u ® u ® w) T = (/""(—x, yi, y 2 , 2/3) 



t - JSS/ 



(3.5) 



where it 



1 
-1 



. This implies that the sign of x does not change the entanglement 



class of p ES (x, yi, y 2 , 2/3). Therefore, it is convenient to restrict ourselves to x > in the 
following. 

IV. SEPARABLE STATES 



In this section we will find a region in the four-dimensional (x, y-\ , yr>.,y?) space, where the 
separable states reside. The calculation procedure is similar to Ref. jl5|. First, we define a 



general form of the fully separable three-qubit pure state by making use of the local-unitary 
transformation, i.e., \ip sep ) = {U\ £g> l/ 2 C§> U%) 1 000) , where 

U 3 =( Aj ~ B n \B^ = 1-\A^. (4.1) 

Second, we map \ip sep ) to the extend GHZ-symmetric state p E (ip sep ) by using a map dis- 
cussed in the previous section. Finally, we maximize x when y 1( y 2 , and y 3 are fixed. 
Combining Eq. (JE3]) and Eq. (|Pjl . p ES (ip se P) is given by Eq. flU]) with 



X = 




l^ilVi- 


\M 2 < 


Ml = 


|A 2 | 2 L4 3 | 2 + (1- 


l^2| 2 )(l- 


\A 3 \ 2 ) 


= 




l^l| 2 )(l- 


\A 3 \ 2 ) 


^3 = 


|Ai| 2 |A 2 | 2 + (l- 


l^l| 2 )(l- 


\A 2 \ 2 ) 



(4.2) 



where \ii = 2yi + ~. In order to apply the Lagrange multiplier method we define 

3 

i=l 

where Aj's are the Lagrange multiplier constants and 0,'s are the constraints given by 



2 • 

3 

x A = x 



6i = 


|a 2 | 2 |a 3 | 2 + (i- 


i^i 2 )(i - 




- Ml 


(4.4) 


e 2 = 




i^i 2 )(i- 


|4>| 2 ) 


- /U 2 




e 3 = 


i^iI 2 ia 2 | 2 + (1 - 


i^ii 2 )(i - 


l^2| 2 ) 


- M3- 





Before proceeding further, it is worthwhile to compare Eq. (14.31) with the corresponding 
equation derived for GHZ-symmetric case at this stage. For GHZ-symmetric case 15j x A for 
the separable states becomes 



x A 



A 1 I|A 2 )IA 3 |V / 1 " l^ilVl " KIM - l^l 2 + A6 (4.5) 
= \A X \ 2 \A 2 \ 2 \A*? + (1 - |A 1 | 2 )(l - |A 2 | 2 )(1 - |A 3 | 2 ) - (V3y + . 

Thus x A in Eq. ( 14. 5 h has a f> Aj symmetry. Thus, the maximum of x occurs when 
I Al| = |^2| = 1-^3 1, which drastically simplifies the calculation. However, as Eq. f 14 . 3 H and 
Eq. ( 14. 4p show, x A in Eq. H4.3[) does not have this symmetry. This is due to the fact that 
the extended GHZ symmetry is less symmetric than the GHZ symmetry. 
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The Lagrange multiplier method generates three equations Jpq- = (i — 1,2,3). Since 
we have three Lagrange multiplier constants, these equations can be used to express Aj in 
terms of \Ai\. Thus, we should determine \Ai\ from only three constraints Oj = 0, which 
yields 

, , |2 V\ ± ^y/Vmyi , . |2 _ y 2 ±2 v /y 1 y 2 y 3 2 y 3 ± 2 y /y 1 y 2 y 3 _* 

\ A l\ - 7, A 2 - ~ ^3 - o ' 

2yi 2 V2 2y3 

Therefore, \Ai\ = \A 2 \ = \A 3 \ doe not hold unless 2/1=2/2 = 2/3- This is due to the less- 
symmetric nature of the extended GHZ symmetry compared to the GHZ symmetry. From 
Eq. (HSj) x max is given by 



1 i (2/1 - 4ym){y2 - tymHys - tym) (A7 , 

Eq. (14. 7p gives a certain boundary in the four- dimensional (x,yi,y 2 ,y 3 ) space, inside of 
which the extended GHZ-symmetric separable states reside. It is worthwhile noting two 
points at the present stage. First, if the term inside the square root in r.h.s. of Eq. ( 14. T}) 
is negative at some point (or region) inside the polyhedron of Fig. 2(a), this means that 
this point is excluded from the boundary. This is similar to — < y < region in the 
GHZ-symmetric case as Fig. 1 exhibits. Second, if the region generated by Eq. (14 ,7p is not 
convex, we should extend it to its convex hull because the set of each entanglement class 
should be convex set. 

Now, we consider several special cases. If 2/1 = 2/2 = 2/3 = y, Eq. (I3.3P gives — i < y < | 
and x max becomes 

Xmax = |(1 - 42/) 3 / 2 . (4.8) 

Since this is not convex (see blue line of Fig. 3(a)), we have to choose a convex hull, which 
is 

Xmax = g ~ 2* (4-9) 

This region is depicted in Fig. 3(a) as a green color. 

As a second example, let us consider a case of 2/1 = 2/2 = — 2/3 = 2/- 111 this case Eq. (13. 3p 
gives — t < y < Then, from Eq. ( 14. 7p x max reduces to 

Xmax = ^(l + 42/) 3/2 . (4.10) 

Since this is not convex (see blue line of Fig. 3(b)) and it is evident that the states with 
x = are separable, we should choose a separable region as a green color region in Fig. 3(b). 



FIG. 3: Separable region for (a) yi = y 2 = y 3 = y, (b) y x = y 2 = -y 3 = y, (c) y\ = y 2 = and 
y 3 = y. For first two cases Eq. (|4.7p generates concave regions (see blue lines). Thus the convex 
hull (red line) for each case is chosen as a separable region. For last case, however, Eq. (|4.7p 
generates a convex separable region (see red line). If y\ = —y 2 = —y 3 = y, the separable region in 
the (x, y) plane becomes the same with (b) upside down. 

If yi = —y 2 = — 2/3 = y, the separable region in the (x, y) plane becomes the same with Fig. 
3(b) upside down. 

As a third example, let us consider a case of (y\, y 2 , y%) = (0, 0, y). In this case Eq. (13. 3p 
gives —7<y< \ and Eq. (14. 7p yields 

X m ax = g(l~%). (4.11) 

The corresponding separable region is plotted in Fig. 3(c) as a green color. Since it is 
convex, we do not need to choose a convex hull in this case. 

Finally, let us consider the positive partial transpose (PPT) states in the extended GHZ- 
symmetric states p ES (x,yi,y 2 ,y3). Taking a partial transposition over the first qubit and 
computing the eigenvalues of the resulting matrix, one can derive the PPT condition for 
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positive x as x < a^, where 

1 2/i + 2/2 ~ 2/3 
a i = T, 7, > 



«3 



1 Vi - 2/2 + 2/3 



Q i 



-2/i + 2/2 + 2/3 



(4.12) 



Performing similar calculation for second and third qubits, it is straightforward to derive a 
PPT condition of p E (x, yi, y 2 , 2/3) as 



x < x max = min(a 2 , « 3 , a 4 ). 



(4.13) 



One can show that the separable regions in Fig. 3 exactly coincide with the region, where 
the PPT condition ( 14+3]) holds. 

V. ROUGH CLASSIFICATION USING THE CLASS-SPECIFIC OPTIMAL WIT- 
NESS OPERATORS 
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FIG. 4: Rough SLOCC classification of the extended GHZ-symmetric states p ES given in Eq. 
(13. 1|) . G, W, and B stand for GHZ class, W class, and bi-separable class. The symbol '+' means 
coexistence. For example, W+G means the coexistence of GHZ and W classes. 
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The Lagrange multiplier method used in the previous section to derive the region for the 
separable states cannot be used to derive the region for the bi-separable states. The reason 
is as follows. If we choose first qubit as a separable qubit such as \i/j btse P) = {G\ <8> G 2 <S> 
G 3 )|0>®^(|00) + |11)), where 



Gj = 1 ' \ , (5-1) 




the mapping from a set of the three-qubit pure states to the set of the extended GHZ- 
symmetric states cannot change the separability of the first qubit because the qubit per- 
mutation is not involved in the extended GHZ symmetry. Since, however, the definition 
of the bi-separable mixed state means a quantum state whose pure-states ensemble can be 
represented as only separable and bi-separable states without restriction to the separable 
qubit, the Lagrange multiplier method used in the previous section cannot be applied for 
deriving the region for the bi-separable extended GHZ-symmetric states. 

In this paper, instead of the Lagrange multiplier method, we use the class-specific optimal 
witness operators 



W bise p\ sep = 1 - 4|GHZ + )(GHZ + |+2|GHZ_)(GHZ_| (5.2) 
1 
2 



>V w \ bisep = ^l-|GHZ + )(GHZ + | 



WghzwM = jl- 2 9 3 — : |GHZ + )(GHZ + |- 3 |GHZ_)(GHZ_|, 
4 Vq - 2v + 4 Vq + 2v + 4 



which are derived in Ref . [U| by using the classification of the GHZ-symmetric states. Here, 
we choose vq = 0.981, which corresponds to the fact that the optimal witness operators yield 
an exact classification to the Werner state 

p w =p|GHZ+)(GHZ+|+(l-p)-l. (5.3) 

8 

The information the class-specific witness Wa\b provides is as follows. Let p be an 
arbitrary three-qubit quantum state. If tr(WA\BP) < 0, this means that p is in A or its 
higher class in the three-qubit hierarchy ScBcWc GHZ. 
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Using Eq. (15. 2p and Eq. (13. ip it is straightforward to show 



tr [WbisepVsep/^ 5 '] = - - (Vl + 2/2 + J/3) ~ 6x 



(5.4) 



tr [W w \ hiS e P p ES ] 
tr [W GU z\wp ES ] 



4 



- (Vi + y2 + 2/3) - 2x 



3(vi + 4) 



2v + 4)(v 2 + 2v + 4) 



A{vl + 4) 



(2/1 + 2/2 + 2/3) 



It is worthwhile noting that Eq. (15 .4p is not dependent on yj (j = 1,2,3) individually, but 
dependent on 2/1 + 2/2 + 2/3- The information we can gain from Eq. (15.41) is as follows. The ex- 
tended GHZ-symmetric separable states should be confined in a polygon (b, f, g, h) in Fig. 4. 
The extended GHZ-symmetric bi-separable states should be confined in a polygon (b, e, g, i). 
The extended GHZ-symmetric W states should be confined in a polygon (a,c,d,g,j). Of 
course, all SLOCC classes should be distributed with obeying the three-qubit hierarchy S 
C B C W C GHZ. This information is pictorially depicted in Fig. 4. The three examples 
discussed in the previous section can be shown to be consistent with this information, i.e., 
all green regions in Fig. 3 are contained in the polygon (b, f, g, h). 

There is another way, which enables us to get a rough classification of the extended GHZ- 
symmetric states p ES . First, we map from p ES in Eq. (13.1 ft to p s in Eq. H2.2I) . which results 



in p (p ). Then, the parameters x and y of p (p ) are 



x 



2 (Pooo.m + Pin,ooo) 
1 

7! 



X 



(5.5) 



Pooaooo + Pin, in 



('2/i + 2/2 + 2/3)- 



Since p s (p ES ) should be lower class than p ES , we conclude 

(i) p ES is a GHZ class if p s (p ES ) is a GHZ class 

(ii) p ES is a GHZ or W class if p s (p ES ) is a W class 

(iii) p ES is a GHZ, W, or B class if p s (p ES ) is a B class. 
This makes a similar (but not exactly same) figure to Fig. 4. 



VI. CONCLUSION 



In this paper we analyze the SLOCC classification of the extended GHZ-symmetric states 
p ES , which is parametrized by four real parameters. The condition for separable states of 

13 



p ES is analytically derived (see Eq. (14.7j) ). The higher classes such as B, W, and GHZ classes 
are roughly classified by making use of the class-specific optimal witnesses and map from 
extended GHZ symmetry to GHZ symmetry (see Eq. (j5.5p ). From this analysis we guess 
that the entanglement classes of p are not dependent on yj (J — 1, 2, 3) individually, but 
dependent on y 1 + y 2 + y 3 collectively. Unfortunately, we do not know how to prove our 
guess from the analytical ground. 

The entanglement classification for the GHZ-symmetric case can be extended to the 
higher-qubit systems. However, analysis of the entanglement classes in the higher-qubit 
systems seems to be much more difficult than that of the three-qubit case. For example, 
the general form of the GHZ-like-symmetric states 3 in four qubit system is parametrized by 
three real parameters in a form 

p s A = (3 [|0000)(1111| + |1111)(0000|] (6.1) 
+diag (an, a 2 , «2, ce 3 , a 2 , a 3 , a 3 , a 2 , a 2 , a 3 , a 3 , a 2 , « 3 , a 2 , a 2 , «i) 

with ol\ + 4a 2 + 3«3 = |. Therefore, total set of the GHZ-like symmetric states should be 
represented by three-dimensional volume in the parameter space. Furthermore, although the 
entanglement classification of the four-qubit pure system is treated in several papers [12|, [17- 
21], their results can be confusing and seemingly contradictory. The worst thing is that the 
entanglement classes of the four-qubit mixed system are not well understood so far and it 
is not clear whether or not they obey the linear hierarchy. We hope to revisit this issue in 
the future. 
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